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SUMMARY 


This report describes the results of a theoretical study of viscous drag 
reduction schemes for potential application to the fuselage of a long-haul sub- 
sonic transport aircraft. The schemes which were examined included tangential 
slot injection on the fuselage and various synergetic combinations of tangen- 
tial slot injection and distributed suction applied to wing and fuselage sur- 
faces. Both passive and mechanical (i.e., utilizing turbo-machinery) systems 

were examined. 

Overall performance of the selected systems was determined at a fixed sub- 
sonic cruise condition corresponding to a flight Mach number of M^ = 0.8 and an 
altitude of 11,000 m. The nominal aircraft to which most of the performance 
data was referenced was a wide-body transport of the Boeing Ihl category. Some 
of the performance results obtained with wing suction are referenced to a Lock- 
heed Star Lifter wing section. 

The results of this study show that very substantial reductions in fuselage 
viscous drag are achievable with tangential slot Injection. However, the drag 
penalties attributable to the components of the baseline design, viz., the turbo 
machinery and ducts, largely offset the reduction in fuselage viscous drag. In 
some cases, a net increase in drag accrues to the system. The drag penalties 
are incurred in the process of reducing the momentum of the captured air used 
for Injection and In pumping this air through the ducts. Perturbations in the 
baseline design, which involved use of ten (10) slots distributed at equally 
spaced intervals over the length of the fuselage, did not Indicate that any 
significant Improvement in the performance of the baseline design was possible 
through minor parametric variations. However, it is clear that if alternate 
designs which avoid the system penalties associated with capturing and pumping 
high weight flow rates of air can be devised, slot Injection, per se, can pro- 
duce significant viscous drag reductions (i.e., as large as 50 ^) with minimal 
complexity and impact on the basic fuselage configuration. 

Alternate designs investigated in the present study involved combinations 
of boundary layer suction on the wing surfaces and Injection on the fuselage, 
and suction and injection combinations applied to the fuselage only. 

-vl 1 i - 



Overall system performance for these designs Is found to be superior to 
the baseline injection scheme, due primarily to reduction of the requisite 
flow rates. The latter are based on theoretical estimates of the suction re- 
quired to maintain a laminar condition in two-dimensional incompressible bound- 
ary layer flow. Whether these theoretical estimates will prove valid under ac- 
tual flight conditions is not known. Nevertheless, with this caveat in mind, it 
is shown that the considered alternate designs offer significant Improvements 
in fuel consumption and/or range characteristics. Fuel load decreases of the 
order of 5 to 17^ and range increases of 8 to 32% are obtained. 



I. INTRODUCTION 


The objective of the present study was to evaluate quantitatively a fuse- 
lage viscous drag reduction system for a representative subsonic aircraft* At 
its inception, this study was structured around a baseline scheme which in- 
volved tangential slot injection through ten (lO) slots at equally spaced in- 
tervals along the fuselage. The study was to include theoretical calculations 
for all system components, including rotating machinery performance and effi- 
ciency and duet losses. 

In the course of this investigation it was found that the baseline scheme 
was incapable of providing the anticipated overall system drag reduction and, 
in some cases, resulted in increased overall system drag. Accordingly, with 
the agreement of NASA, alternate schemes were examined. These involved the • 
use of boundary layer suction on both fuselage and wing surfaces in various 
combinations with slot injection. These schemes have been found to provide 
significant improvement in overall performance. 

The efficacy of boundary layer suction In reducing viscous drag resides 
in the generally recognized principle that suction can stabilize a laminar 
boundary layer. As a result, the skin friction on a surface through which 
suction is applied can be reduced to a small fraction of its value for the 
naturally turbulent boundary layer on the same surface. 

Experimental evidence exists demonstrating the abl 1 i ty of suction to main- 
taining a laminar flow for the conditions of Interest in the present investi- 
gation under ideal conditions. For the subsonic case the results . of Refer-, 
ences (1) and (2) may be cited. The results of Pfenninger (Reference 3) indi- 
cate that similar results can be achieved for supersonic flow for a length 
Reynolds number up to 5 x 10 . The feas i b i 1 i ty of suppressing separation and 
maintaining laminar flow In and downstream of interactions with weak incident 
shock waves has been demonstrated by the experiments of Groth et al (Reference 
. On the other hand, experience shows that there are formidable problems 
associated wi th utilizing laminar flow control when nonuniformities in the suc- 
tion distribution and/or surface geometry are present. As will be seen later, 
the schemes which are examined here inherently involve such nonuniformities. 



In addition, the fabrication, operation and maintenance problems associated 
with a suitable porous surface have not been addressed in this study. Ac- 
cordingly, while the results presented here delineate the potential of these 
schemes, they should be considered provisional in terms of application. 

In the present study the potential benefit in performance due to suction 
techniques was examined in two distinct ways. In the first of these, suction 
was applied to essentially constant pressure surfaces (fuselage) to stabilize 
a laminar boundary layer so as to prevent transition and the associated in- 
creases in viscous shear stress on the surface. In the second approach, sue-*' 
tion was applied to wing surfaces experiencing adverse pressure gradients 
with the aim of preventing separation. In these circumstances, improvement 
in aircraft aerodynamic performance would accrue both from reduction in vis 
cous shear as well as from improved L/D characteristics of the wing. 

Theoretical estimates of viscous drag reduction due to tangential slot 
injection were provided by the NASA Langley Research Center in accordance 
with contractual agreement. Corresponding estimates of laminar boundary 
layer behavior with suction were generated by ATL employing various approxi- 
mate schemes which are described in subsequent sections of this report. 

The overall performance of the selected systems was determined at a fixed 
subsonic cruise condi tion corresponding to a nominal flight Mach number of 
M = 0.8 and an altitude of 11,000 m. The aircraft configuration to which 

oa 

most of the performance data was referenced was the Boeing 7^7* Some of the 
performance results obtained with wing suction are referenced to a Lockheed 
C-l4l Star Lifter wing section. In this connection it must be emphasized 
that no optimization in terms of aircraft configuration was attempted in the 
present study. Accordingly, the results obtained can probably be improved 
by sppropriate changes in configuration. Recommendations in this regard are 
presented in the last section of this report. 
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II. LIST OF SYMBOLS 


a 

A 

b 

B 


c 



L 

L. 

I 


L/D 

(L/D) 


(L/D)/(L/D) 
flow area 
slot width 
AW/W 

o 

chord length 

local skin friction coefficient with suction or injection 

average skin friction coefficient with suction or injection 

local baseline skin friction coefficient 

average baseline skin friction coefficient 

pressure coefficient 

suction flow coefficient - p v /p V 

s s 

fuselage diameter 

net drag with suction or injection 

friction drag on surface affected by suction or injection 

baseline total drag 

baseline fuselage drag 

basel i ne fuel load 

fuel load with drag reduction 

slot height 

turbo-machine drag per unit flow rate (see page 37) 

line loss recovery factor 

length of suction or injection interval 

streamwise distance to first slot 

overall length of fuselage 

baseline maximum lift-to-drag ratio 

maximum 1 i ft-to-drag ratio wi th drag reduction 



M 

N 

P 

q 


r 


r 

R 

I 

R 


Re 


f 


V 

s 

V 


w 


W 


(U 

e 

bl 

o 

X 

X 


X 


y 


Y 
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n 


c 


e 

X 


Mach number 

mass flow function (see page 66) 

pressure 

dynamic pressure 

fuselage radius 

r/c 

baseline range 

range with drag reduction 

Reynolds number based on fuselage length and free stream conditions 
suction velocity 
axial velocity 

aircraft weight increment per unit surface area 

mass flow rate 

W - F 
o 

gross weight of ai rcraft 
streamwise coordinate 

skin friction reduction parameter (see page k3>) 

normal coordinate 

specific heat ratio 

boundary layer thickness 

displacement thickness 

compressor efficiency 

turbine efficiency 

momentum thickness 

slot injection parameter ~ P;V./p 

J J 


M 


viscosity 


K suction parameter - C p V y/p 

p density 


Subscripts 

c compressor Conditions 

j slot injection conditions 

s suction conditions 

1 turbine entrance conditions 

2 turbine exit conditions 

3 compressor exit conditions 

^ compressor entrance conditions 

“ free stream conditions 


Superscri pts 
( ) average values 


sonic conditions 



III. PRELIMINARY CONSIDERATIONS 

A. Baseline Aircraft Drag - This study was conducted for typical CTOL 
cruise flight conditions and a fuselage shape representative of current 
long-haul subsonic transports. Most of the results of this study are pre- 
sented in terms of net drag reduction for the various schemes as a percent 
of a reference total drag D^ corresponding to the selected basel ine ai rcraf t. 
For the present purpose 0^ has been taken to correspond to the total drag of 
a wide body transport (viz., a Boeing 7^7) cruising at an altitude of 11,000 m 
and M^ = 0.82. In order to estimate this parameter the following approxima-' 
tions were employed. 

For the cited flight conditions the unit free stream Reynolds number is 
approximately 6.2 million/m so that, based on a fuselage length of approxi- 
mately 67 meters, a fuselage Reynolds number. Rep, on the order of M x 10 
prevails. An average skin frictior. coefficient based on this Reynolds number 
can be obtained from the Prandtl-Schl Ichting correlation (Reference 6) yield- 

JUJL 

ing"'* 

Cp = 0.455 (log Rep)"^*^® = .00175 

Taking the diameter of the fuselage to be 6.7 meters the total wetted area 

2 

is approximately 1400 m . Accordingly, Dp , the fuselage drag in the absence 
of any drag reduction effects is approximately 27,000 newtons. 

An average skin friction coefficient for the wing is estimated to be 
given by the above skin friction law but for a length Reynolds number based 


*Characteristic dimensions of the various aircraft considered here are taken 
from Reference (5). 

**Compressibi 1 i ty effects would reduce this value by approximately 10% for 
the adiabatic wat 1 case. However, this value is considered sufficiently 
accurate for the present purpose. 
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on an average chord length of 90 meters. This value (.00231) is applied to 
an exposed surface area of around 81»0 m . The skin friction drag of the wing 
is, therefore, 21,000 newtons. Adding 10^ to adcount for skin friction on the 
empennage and engine pods and profile drag the total zero lift drag of the 
baseline aircraft is estimated to be 53,000 newtons. 

The airplane cruises at maximum (L/D). The drag corresponding to this 
condition is twice the zero lift drag. Thus, the cruise drag of the baseline 
aircraft, D^, is estimated to be 106,000 newtons. 

B. Fuselage Geometry and Pressure Distribution - As indicated earlier, 
calculations of the turbulent shear distribution on the fuselage in the pres- 
ence of tangential slot injection were carried out at the NASA Langley Re- 
search Center. The numerical finite difference technique due to Beckwith and 
Bushnell (Reference 7) was employed for this purpose. 

In order to implement this methodology a body geometry and corresponding 
pressure distribution were needed. The fuselage was approximated by the quasi- 
ellipsoid of revolution depicted in Figure (l). In the absence of experimental 
data an estimate of the pressure distribution was made using the method of 
Reference (8). This is depicted by the solid line shown in the lower portion 
of Figure (1). Note that this distribution indicated the existence of sub- 
ambient static pressures over much of the fuselage surface as well as singu- 
larities at the fore and aft stagnation points. 

To simplify the numerical procedures, the alternate pressure distribution 
shown by the dashed line was substi tuted wi th the agreement of the contract . 
technical monitor. 

C. Wing Geometry and Pressure Distribution - Accurate calculation of 
the potential benefit of wing suction requires both three-dimensional tran- 
sonic inviscid analysis and three-dimensional laminar boundary layer analysis. 
The latter should include the s imul taneous effects of suction and rapid 
streamwise variations in pressure as well as, possibly, stabil i ty cons 1 dera- 
tions. Clearly, such a detai led approach would be beyond the scope of the 
present design study. 
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To provide at least a rough estimate of the effect of suction on the 
baseline configuration a two-dimensional approach was adopted. Within this 
Idealized framework the essential feature of a super-critical airfoil with 
the attendant weak shock was retained. In the absence of any direct informa 
tion on the Ihl airfoil the inviscid pressure distribution over a C-Ul air- 
foil section, as obtained from Reference (9), was employed. This pressure 
i but i on for both upper and lower surfaces is shown in Figure (2). It 
corresponds to flow at = 0.76 and an angle of attack of 0.95°. These 
values are not precisely those of the baseline configuration, particularly 
the angle of attack. Nevertheless, the resulting pressure variation should 
be representative in terms of the presence of a normal shock on the upper 
surface. 


The airfoil section geometry is shown in Figure (3). 







FIGURE 3. WING GEOMETRY 





IV. BOUNDARY LAYER CALCULATIONS 


A. Constant Pressure Laminar Results - Skin friction reduction due to 
uniform suction at constant pressure was estimated directly from the incom- 
pressible results given by Schlichting (Reference 6). These results are 
summarized in Figure (A). For the present application the parameter of 
interest is the net reduction in average skin friction coefficient, Cp, re- 
lative to the baseline turbulent value, Cp . The dependence of this ratio on 
the relative suction rate C^ has been deduced from the results shown in Fig- 
ure (4) and are presented in Figure (5) with the length Reynolds number as a* 
parameter. 

Indicated in Figure (5) is a "cutoff" value of C^ = .00012. This value 
corresponds to the minimum needed to insure the maintenance of laminar flow 
as established by the stability considerations outlined in Chapter XVII of 
Reference (6). - 

Finally, we note that these incompressible results should be reasonably 
accurate for the present high speed application since only the ratio of skin 
friction levels is involved. 

B. Laminar Results with Pressure Gradient - The analytical method uti- 
lized for these calculations is described in detail in Appendix A. It em- 
ploys the momentum integral technique due to Torda (Reference 10) in conjunc- 
tion with Thwaites method (Reference IT) to permit initiation of the calcula- 
tion at a stagnation poi nt. 

The numerical computation scheme based on this method can be exercised in 
two distinct ways. The more general option accepts arbitrary distributions 
of suction (or blowing) and pressure gradient and determines all boundary layer 
characteristics. The alternate option imposes the condition that the boundary 
layer thickness is constant. In this case, for an arbitrary variation of ex- 
ternal pressure, the computation yields the requisite suction distribution as 
well as the corresponding variation of all other boundary layer properties. 

The latter option is useful in terms of providing a mechanism by which, in an 
approximate sense, the preservation of a laminar flow can be assured. 


12 





3 



FIGURE 5. REDUCED SKIN FRICTION COEFFICIENT WITH SUCTION 
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In applying this analysis to the baseline wi ng configuration the follow- 
ing objectives prevailed; the suction distribution should insure that the 
state of the boundary layer remain laminar; the suction distribution should 
prevent separation under the influence of the inviscid pressure distribu- 
tion associated with the baseline configuration; the suction distribution 
should yield the minimum skin friction consistent with the previous require- 
ments. Toward this end the suction distributions on the wing were determined 
as follows. 

In regions of constant pressure and in regions of favorable pressure 

gradient (which tends to suppress the rate of growth of the boundary laye*-) , 

the relative suction rate was maintained constant at the optimum value 

C = .00012. In principle, this can be expected to preserve the laminar 
s . . 

boundary layer state. In regions of adverse pressure gradient the variation 

of C was computed by requiring that the boundary layer thickness remain con- 
s 

Stant at the value associated with the start of the pressure rise. The cor- 
responding momentum thickness is found to decrease through this region so 
that, here again, it can be anticipated that the laminar state will prevail. 

The results of these calculations are shown in Figures (6) and (7). As 
can be seen, very modest increases of suction over the optimum value are re- 
quired to prevent separation on the lower surface of the wing. On the upper 
surface, of course, a very large "spike" of auction intensity is needed at 
the ^Q% chord station and at the trailing edge to maintain an attached laminar 
flow. Note, however, that although the maximum value of is on the order of 
kO times the optimum rate it is still less than ,1^ of the unit free stream flow 
rate. 

The corresponding distribution of skin friction is shown in Figure (7) • 
these laminar distributions have been compared with turbulent estimates made 
using the method described in Reference (12). Note that for the upper surface, 
separation is predicted at the ~J0% station and that the turbulent shear is 
assumed to be vanishingly small thereafter. 

The pertinent results needed for the performance calculations are the net 
reduction In average shear for the enti re wing and the average suction rate. 
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These data have been computed from the above results and are summarized be- 
low In Table I. The results indicate that the wall shear is reduced to ap- 
proximately 25% of the turbulent values by application of suction at a rate 
on the order of twice the optimum value for a flat plate. 


TABLE I 

SUMMARY OF VISCOUS DRAG REDUCTION RESULTS ON WING 



Upper 

Lower 

Total 


Surface 

Surface 

Wing 

Average turbulent shear Cp 

i 

.00175 

.00283 

.00229 

Average laminar shear Cp 

.000575 

.000375 

.000475 


.343 

.133 

.238 

Average suction parameter C^ 

.000375 

.000185 

.00028 


C. Turbulent Results With Slot Injection - Turbulent boundary layer so- 

Ijutions with tangential slot injection were obtained by F. G. Howard at the 

i 

NASA Langley Research Center . These were carried out for the fuselage con- 
figurations previously shown in Figure (l). The first slot is located at 
x/L^ = 1/11 and subsequent slots are placed at intervals of Ax/L^ - 1/11, 
to a maximum of 10 slots. 

The assumed surface pressure distribution up to the first slot is shown 
in Figure (8); from the first slot to the end of the fuselage the pressure co 
ejfficient was assumed constant and equal to zero. The numerical method of 
Reference (1A) was used to calculate the boundary layer characteristics up to 
the first slot, assuming a fully developed turbulent boundary layer from the 


*A more complete presentation of these results may be found in Reference (13) 
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nose. The boundary layer thickness just upstream of the first slot is 
5 = 3*6 cm, the displacement thickness is 6 = .57 cm, and the momentum 

th icKness is 9 = .35 cm. The boundary layer velocity profile was then com- 
bined with estimated slot exit velocity profiles having a shape similar to 
those measured in Reference (15), and an average Mach number of = 0.2. 

The resultant complete velocity profile was used as input for the slot in- 
jection code of Reference (7). The resultant velocity profile at the first 
slot (for h = 7.62 cm) is shown in Figure (8). The slot to free stream to- 
tal temperature ratio was assumed constant and equal to 0.9895. The varia- 
tion of skin friction and velocity profiles downstream of asingle slot is - 
indicated in Figure (9). 

The numerical finite-difference solution of Reference (7) was modified 
so that the effect of multiple slot injection on the fuselage skin friction 
could be determined. The local skin friction coefficients (C^) obtained 
downstream of one, three, five and ten slots (slot height (h) of 7.62 cm) 
are compared with the local Skin friction coefficient on the fuselage with- 
out slots in Figure (10). The local skin friction reduction with only one 
slot is significant when compared to that without slot injection. The bene- 
ficial effect of the slot injection is most pronounced immediately downstream 
of the slot exit and diminishes with increasing distance downstream from the 
slot; this occurs because in the near slot region the wall friction is influ- 
enced only by the slot flow while further downstream mixing between the high 
momentum boundary layer flow and the relatively low momentum slot flow in- 
creases the wall shear. The effect of slot height on skin friction varia- 
tions with downstream distance is demonstrated in Figure (II) for the ten slot 
case. For these calculations the slot velocity profiles were scaled by the 
slot height so that the slot mass flow varied in proportion to slot height. 

As in the case of the laminar results, the net reduction in average skin 

friction relative to the baseline value, C^-ZC^ , is needed for the performance 

• . . ^ 

calculations. For this purpose the average wall shear of the fuselage was 

estimated from the data shown in Figures (10) and (11) and normalized with 
respect to the corresponding average for the no slot case. The resulting 
variation with number of slots and slot height is presented in Figure (12). 
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FIGURE 9. VELOCITY DEVELOPMENT AND SKIN FRICTION BEHAVIOR DOWNSTREAM OF A 
SINGLE SLOT 
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The improvement with additional siots and increasing slot height, both of 
which correspond to increasing mass flow, is apparent. It is clear from 
Figure (12) that large reductions (s. 50« in viscous drag are available 
through the use of slot injection systems. 

The results shown in Figures (10), (11) and (12). suggest that the skin 

friction reduction is' improved by increasing the number of injection 
slots but at a diminishing rate (for constant slot spacing). One probable 
reason for this is that slot location is very important; for the present 

study, the most forward slot is the most effective and the .most rearward^ . 

slot is the least effective. Two advantages of a forward slot location m 
the present study are (l) the no injection skin friction level is high and 
(2) the boundary layer is thin; slot effectiveness for local skin friction 
reduction is improved at low ratios of boundary layer thickness to slot 
height (References 16 and 17) • Forward slot injection offers an obvious 
additional advantage in that the drag reduction occurs over a larger area 
of the aircraft. This effect of slot location is illustrated when the 
following comparison is made from Figure (12); consider the case of ten 
slots with h = 3.81 cm compared with the case of five slots with h = 7.62 cm 
(t(ie first slot in each case is located at the same position). Although the 

total mass flow from the five larger slots would be the same as that for the 

ten 3.81 cm slots, the skin friction reduction is H% greater for the five 
slot configuration than for the ten slots (see comparison in Figure 12). 

The effect of injection Mach number has also been investigated. The 
results are summarized in Figure (13), which includes variations in slot 
height as well as injection Mach number. It appears that reduction of the 
injection Mach number can provide further reduction in the skin friction 
fcir the same total mass flow (i.e., compare points A and B in Figure 13). 
However, as wi 1 1 be seen subsequently, substantial drag penalties are in- 
curred- in the process of capturing the required mass flow and reducing its 
average Mach number to the level desired at the slot locations. The total 
mass flow injected through all the slots is captured downstream of the last 
slot, in the proposed baseline design, and reprocessed. However, due to 
the mixing which- occurs in the fuselage boundary layer downstream of each 
slot, the captured flow is more energetic than desired. The excess energy 
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FIGURE 13. SKIN FRICTION REDUCTION EFFECTIVENESS AS A FUNCTION 
OF SLOT HEIGHT AND INJECTION MACH NUMBER 







is removed in the turbo-machinery used for pumping the Injected flow, which 
will be discussed subsequently. This problem Is aggravated by reduction of 
the slot height, which increases the mixing and energization processes, as 
indicated by the velocity profiles at a station near the end of the fuselage 
shown in Figure (l4), and by reduction of the injection Mach number. 






. V. DRAG REDUCTION RESULTS 

A. Fuselage Slot Injection - The scherne for utilrzing the drag reduc- 
tion potential of slot injection on the fuselage in conjunction with a 
turbine/compressor illustrated in Figure (15). The turbine processes the 
boundary layer air at the end of the series of slots and returns the air- 
through an annular duct around the fuselage from which it is injected to the 
slots. Since the flow is injected at a low velocity in order to reduce skin 
friction, there Is a drag associated with the turbine flow. The power gen- 
erated by the turbine is absorbed by the compressor which also processes 
boundary layer air. The compressor air is discharged at a high velocity pro- 
ducing thrust which partly offsets the drag of the turbine flow. 

As noted in Figure (I 5 ) there are two ways in which the boundary layer 
air can be processed. In the first way, the compressor can handle the "inner" 
or lower velocity flow near the surface whi le the turbine handles the "outer" 
or high velocity flow. Alternately, the location of turbine and compressor 
are reversed, with the turbine handling the"inner" flow and the compressor 
handling the "outer" flow. The arrangement which would be selected is the 
one giving the smallest net drag when processing the boundary layer air pro- 
duced by the injection. 

From the velocity profile data such as those presented in Figure (1A) , 
the average velocity as a function of boundary layer flow can be determined. 
Here, the average velocity means the average "momentum'* ve loci ty defined by 

V = n -^vdw 

w . 

whiich is judged to be the most pertinent average for use in calculating the 
dr|ag of the turbo-machine flows. For convenience of application, the average 
velocity is determined as a function of W/Wj , the boundary layer flow with 
respect to the injected flow, being the total flow injected from all the 
slots. The results are presented in Figures (l6a) and (l6b). 

The average velocity enters into the drag calculation of the turbo- 
machines in a straightforward manner. The drag of the turbine flow is given 
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(b) Turbine Handling Inner Flow 


FIGURE 15. SCHEMATIC OF SLOT INJECTION SYSTEM 











by 


turbine drag - ^ " ^2^ 


where V is the average velocity in the boundary layer flow handled by the 
turbine and V is the discharge velocity or the slot injection velocity. 

This formulation of drag implies that the pressures at intake and discharge 
are the same (equal to free stream pressure), a basic premise of the boundary 
layer calculations which is very nearly true in the real case. 

In reducing the velocity from V, to l/^. the turbine performs work Which 
is transmitted to the compressor. The compressor utilizes the work to in- 
crease the velocity of the flow it handles from the intake value to the 
discharge value Vj. In this case, the compressor flow results in a thrust 

g i ven by 


compressor thrust = (V^ “ 

where is the average velocity of the boundary layer flow handled by the 
compressor. 

The combined turbo-machine drag is then given by 
turbo"machi ne drag = W^I 


where 


W 


-V, - (V3- V^) 


Then , for a g i ven 


In application, the quanti ties V^, and are known, 
value of W /W , the value of can be found by means of a calculation which 
equates the pLer output o^the turbine to the power input of the compressor 
and then transforms the power input to an increased energy of the flow which 
is converted to an Increased discharge velocity V^. 

In formulating the production and conversion of power, conventional tur- 
bine and compressor efficiencies are included. In addition, ducting losses 
are taken into account by means of a loss factor (recovery factor or ratio of 





total pressures) which accounts for the entire loss in the ducting, both 
ahead of and behind the machines. The loss factor is especially important 
in the case of the turbine because of the turning of the flow and the long 
duct length upstream for injection through the slots. 

To convert the average velocity data of Figure (l6) to values of aver- 
age velocity for the turbine and compressor, it is necessary to fix the re- 
spective flows. For the turbine, the flow is fixed a priori since it is 
equal to the injected flow. However, the compressor flow may be chosen 
arbitrarily with the final choice deferred until the effect. on jet drag re- 
duction is determined. In the case where the compressor flow is the "inner" 

flow (see Figure 15) the average velocity V, for a given value of W /W is 
- T c t 

the value of V read at W/V/. = W /W . Then, since the turbine flow is the 

J C c 

same as W., the turbine flow extends from W/W. = W /W to W/W = 1 + W /W . 
TU , J c t j c t 

The average velocity for this portion of the boundary layer flow is given by 


. = (1 + V/yw^) - (W^/W^) V.2 

where is the value of V at W/W = 1 + W^/W^ and .V^ is the value of V at 

W/W. = W /W . c 

J c t 

The average velocities for this case are shown in Figures (17a) and (17b) 

as a function of W /V/ . 

c t 


In the case where the turbine flow is the "inner" flow, the average ve- 
locity is determined as the value of V at W/W. = 1. For a given value of 
W^/W^, the compressor flow extends from W/W. = 1 to W/W. = 1 + W /W and its 
average velocity is given by 



0 


* h 


w /w, 

C t 



where V is the value of V at W/W. = 1 and V., is the value of V at W/W. = I 
' J 2 j 

+ W^/W^. The average veloci ties for this case are shown in Figures (l8a) and 

(18b). 
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Consider now an airplane with a baseline total drag and fuselage 
friction drag D_ . With injection, the fuselage drag becomes 


Dp = X Dp 


where X is a factor representing the effect of injection on skin friction. 

The net drag of the airplane with the injection system, including the 
reduction in skin friction, the drag due to the turbine and the thrust due 
to the compressor, i > given by. 

D = D - (1 - X) + W.I 

o J ' 


The injection flow produces a reduced skin friction drag, defined by the 
factor Cp/Cp , over the portion of the fuselage influenced by the slots. 

The value of'x is not exactly the same as Cp/Cp^ because injection starts 
at some distance from the nose of the fuselage so that there is a small por- 
tion of the fuselage where the drag remains unchanged. In this region assume 
that the friction coefficient varies inversely with distance to the 1/5 power 
and that the wetted area is proportional to the length. Then the value of X 
is given by 


X = Cp/Cp + (1 - Cp/Cp ) 

where 

L. distance to first slot 

i 

L length of fuselage. 


Application of these results to the basel i ne conf i guration is made assum; 
ing the slots to extend completely around the periphery of the fuselage. For 
the cases which were examined the flow rates Wj are 
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V. 5 


The results, which are presented as curves of (net drag with in- 
jection/drag without injection) vs. (compressor flow/turbine or in- 

jection flow) for various values of (line loss factor of injection flow), 
are based on turbine and compressor efficiencies equal to 0.9. This would f 

tend to give results which are slightly optimistic. In the same spirit, 
the line loss factor of the compressor flow was taken to be unity. 

V 

Results for the largest slot height considered (h = cm) are given I 

in Figures (I9a) and (19b) for the compressor handling the "inner" and "outer" 
flows, respectively. Note that the figure includes a line representing the 
effect of injection on skin friction only, given by ^ 

D = D - (1 - X) Dp 

In addition, a line representing results without the turbine/compressor ef- | 

feet Is included. This is given by 

i 

D = D - (1 - X) Dp + W. (V. - V.) 
o J I Z 
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FIGURE 19a. NET DRAG WITH SLOT INJECTION, COMPRESSOR HANDLING INNER 


FLOW; h = 15.2^1 cnij Mj = 0.2 
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FIGURE 19b. NET DRAG WITH SLOT INJECTION, TURBINE HANDLING INNER 


FLOW; h = 15.24 cm, = 0.2 





It is apparent that the turbine/compressor is quite effective in reducing 
net drag although not to an extent sufficient to produce a positive result 
(D/D^ less than unity). Comparison of the results given in Figures (I9a) 
and ( 19 b) also indicates that slightly better results can be achieved with 
the turbine processing the lower momentum "inner" flow and that in neither 
case is the compressor flow of critical importance. These trends prevailed 
for all of the configurations examined. 

Additional results for the smaller slot heights are presented in Fig- 
ures ( 20 ) and { 21 ). Furthermore, to demonstrate the basic soundness of the . 
concept of having the turbo-machines process boundary layer flow, some cal- 
culations were made with the machines handling free stream flow (V^ “ “ '^00 

The comparison shovm in Figure (22) demonstrates the potential benefits. 


Comparison of the results shown in Figures (20) and (21) with those pre' 
sented in Figure (19) indicates that the smallest slot height, despite its 
inferiority in terms of reducing viscous drag, gives the best net perfor- 
mance. This is, of course, a reflection of the large penalty in momentum 
drag which increases directly with total injector flow rate and, therefore, 
with slot height and/or injection Mach number. It should be noted also that 
the effect of the loss factor is greatest for the largest slot height and 
decreases systematically with decreasing height. This is, again, an indica- 
tion of the large penalty In drag associated with increased flow rates. 


As can be seen in Figure (21), even for a loss-free system (k^ = l) , 
net drag reduction is not attainable with Mj = 0.2. Accordingly, additional 
calculations at a lower jet Mach number (Mj = O.I) were carried out for the 
small and intermediate slot heights. These results are presented in Figures 
( 23 ) and (24). Again, even with the decreased effectiveness of slot injec- 
tion in reducing wall shear (c.f. Figure 13) some improvement in overall 
performance is attained due to the lower mass flow rates involved. 


Although the results shown in Figures (23) and (24) indicate some net 
drag reduction the actual values are quite small amounting to only 1 or 2^ 
for realistic values of line loss factors. Further reduction in jet Mach 









































number, or injector flow rate ? 
cant improvement since the wall 
feet is demonstrated in Figure 
goes through unity at Mj -> 0. 


n general, is not likely to provide signifi- 

shear ratio would approach unity. This ef- 

( 25 ) in which the curve of D/D versus M. 

o J 


In summary, these results indicate that the baseline design, which makes 
exclusive use of multiple slot injection and captures the required mass for 
injection downstream of the last slot, will provide, at best, marginal drag 
reduction. It appears that alternate schemes are needed, possibly combining 
slot injection with other drag reduction methods, to provide significant 
improvement in performance. Several alternate schemes are examined in the 
subsequent sections; however the considered schemes do not encompass all pos- 
sibllties by any means. 

B. Combined Fuselage Suction and Slot Injection - This scheme consists 
of a passive system requiring no pumping or turbo-machinery. The system 
depends for its operation on the pressure difference which exists between 
the upper and lower surfaces of the fuselage at a small angle of attack. The 
high pressure lower surface is composed of a porous or slotted surface through 
w^ich suction takes place. The flow is ducted to the top surface of the fuse- 
lage where it is discharged by slot injection. The physical arrangement 
would consist of a series of 10 suction surfaces encompassing the lower half 
of the. fuselage, each of 6.1 m length. At the end of each section, the sue- 
,tion flow is collected and discharged through a slot on the top of the fuse- 
lage, thus comprising a system of 10 slots at 6.1 m intervals. 


As indicated in Figure (5) suction at an excessive rate can produce an 
ihcrease in skin friction rather than a decrease. The largest drag reduc- 
tion is obtained with the "optimum rate" corresponding to the stability 
limit for sustaining a laminar boundary layer. 

The injection at low velocity on the top of the fuselage produces a re- 
duction in skin friction. With the flow fixed by the suction requirement, 
the question arises about the best way to distribute the flow on the top of 





the fuselage. The same flow may be injected with a larger slot height by 
limiting the peripheral extent of the slot. The skin friction reduction 
increases as slot height increases, but the extent of surface affected by 
injection decreases. The best slot height is the one which gives the 
greatest net gain from the two opposing effects. 

The following simple analysis shows that the best slot height is the 
smallest slot height; that is, the flow should be distributed over the 
largest possible lateral extent of surface area. For a given injection 
flow, W, the slot height h and the slot lateral extent b are related by 

W = p . V . hb 
J J 


For the length of the slot coverage L, the area of surface affected by in 
jection is bL. The friction drag on this surface is given by 


Dp = (Cp/Cp ) Cp q. bL 

i i 

where Cp/Cp is the reduced skin friction factor due to injection. 

r r * , 


The drag reduction is then given by 


ADp = (1 - Cp/Cp ) Cp^ 


bL 


I I 


which may also be written 


ADp = 


(1 - Cp/Cp ) Cp LW 
i i 

p. V. h 
J J 


Since C,. , q , L, W, p. and V. are all constant, the drag reduction is written 
r . “ J J 


ADp = constant 


(1 - Cp/Cp ) 

^ i 


It is apparent that the maximum drag reduction is obtained when the value of 
(1 - Cp/Cp )/h is a maximum. The data in Figure (13) are converted to this 
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expression with results presented in Figure (26) showing that the largest 
drag reduction corresponds to the smallest slot height. 

The suction surfaces are arranged in 10 sections of 6.1 m length on the 
bottom half of the fuselage. At the end of each section, the suction flow 
is collected and ducted around the fuselage along which it is injected in a 
slot encompassing the top half of the fuselage. Thus, the slot injection 
consists of 10 slots at 6,1 m intervals. This arrangement does not present 
serious difficulties with regard to internal ducting, especially in view of 
the low flow rates characteristic of suction requirements. It is estimated 
thiat with a 1/2 inch gap around the bottom of the fuselage, the suction flow 
of: each 6.1 m section can be handled with a pressure drop of one-half of 
one percent. Including the ducting loss and the required velocity head for 
injecting the flow through the slot on the upper half of the fuselage, the 
tptal pressure drop will be around 4 or 5 percent. Such a pressure differen- 
tial between the lower and upper part of the fuselage can be obtained with 
a small angle of attack. 


The flow through the suction surface is 


W ^ C p V ^ L 

S oo 00 2 


and the flow through the injection slot is 


W = X p V ^ h 

CO OO 2 


Hence, the slot height is given by 


h = CgU 


The value of determines the skin friction reduction factor due to suction 
) and the value of h determines the skin friction reduction factor 

F F • 

:1s . 

due to injection (Cp/Cp ) . It should be noted that in order to determine 

‘ j V 

(Cp/Cp ) , it is necessary to fix a Reynolds number. The value used in the 


FIGURE 26. PARAMETER DESCRIBING DRAG REDUCTION AS A 
FUNCTION OF SLOT HEIGHT 
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present application is 10 , roughly corresponding to the average fuselage 
length. 

The drag reduction due to suction is given by 


= (1 - (Cp/Cp ) ) Cp q 


Trd 

oo 2 


and that due to injection is 


ADj = 0 - (Cp/C ) ) Cp^ 


the drag increase due to the momentum change of the flow is given by 

= W (V^ - V.) 
m “ J 

which in terms of the suction flow coefficient is written 

AD = 2C q ^ L (1 - V./Vj 
•m s 2 J " 

If it is assumed that the injection density is the same as the free stream 
density (which is very nearly true) then the velocity ratio is the same as 
the injection flow parameter 


M./M = X 

J 00 

The net drag reduction is given by 


AD = AD + AD. - AD 

s J m . 

The original skin friction drag without suction or injection, determined by 
both the upper and lower halves of the fuselage is given by 


D_ = Cc- q TTdL 
F. F, 

I I 


Hence, 


i - 


AD 

D. 


(Cp/C, ) 1 - (Cp/Cp ) 

_ Ls_ + — LL 


(1 - X) 


'F. 


The values of A/D^ determined as a function of the suction flow co- 

efficient C are shown in Figure (27), calculated with Cp = ,00175 and, 

^ i 

X r= .24^. The figure shows the contribution of each of the three compo- 
nents of drag, indicating that suction contributes the largest part of the 
drag reduction. 


In order to evaluate the effect of drag reduction on airplane perfor- 
mance, it is necessary to establish the drag reduction with respect to the 
overall drag. In the application treated earlier, the overall drag was 
established as equal to four times the fuselage drag or =?= A x Dp . This 
relationship will also be used here. However, in order to do so, it is 

necessary to establish the realtionship between Dp and Dp . 

o i 


As in the earlier application, it is assumed that there is an Initial 
fuselage length of 6.1 m which is not treated by suction or injection. The 
drag Dp excludes this length since it represents the original drag on that 
portion of the fuselage treated with suction or injection. Adopting the 
same assumptions used before (turbulent skin friction coefficient varying 
inversely with distance to a .2 power and surface area varying directly with 
distance), the relationship between Dp and Dp is determined as 


F. L. 

= 1 - (^) 

F o 

o 


.8 


With L. =6.1 m and L =67 m, 
I o 


F. 

I 


= .853 


The net drag is found from 


*This corresponds to slot injection at Mj =0.2. All performance estimates 
for this scheme are made using (Cp/Cp_) values corresponding to Mj = 0.2 
since this provides greater reduction in wall shear than the Mj =0.1 case. 










1 


= 1 - .213 




with results as shown in Figure (27). 



It is apparent that the greatest net drag reduction is obtained with 
the minimum suction flow. The injection slot height corresponding to this 
condition is extremely small, approximately one-tenth of an Inch. It is 
possible that a practical installation would require a larger slot height, 
obtainable by restricting the lateral extent of the slot. In this case, 
the contribution to drag reduction would be less than for the fully extend- 
ed slot. However, the slot injection contributes a minor part of the net 
drag reduction so that a larger slot height would not seriously affect the 
net reduction. 

Figure (27) shows that the minimum reduced drag ratio corresponding to 
the minimum suction flow is around .91 {S% drag reduction). The interpreta- 
tion of this ratio in terms of mission performance (range Increase with the 
same fuel load or fuel load decrease with the same range) is presented in a 
later sect ion . 

C. Fuselage Suction - The scheme of combined suction and injection in- 
dicated that the contribution of injection to the net drag reduction is very 
small. The real advantage of the injection is that it takes place in a low 
pressure region on the top half of the fuselage, thus providing a means of 
discharging the suction flow without tte need of pumping or turbo-machines. 

The disadvantage of the scheme is that it precludes achieving the larger drag 
reduction which would result if the upper part of the fuselage were also 
treated with suction. The use of suction over both halves of the fuselage 
would nearly double the reduction in skin friction and would more than off- 
set the drag due to the pumping required. 

, The use of turbo-machines to perform the pumping is neatly adaptable to 
the scheme of suction over the entire fuselage perimeter, permitting a system 
which eliminates the need for complicated ducting. As illustrated in Figure 


59 



(i8) the suction surface forms a shell around the pressurized fuselage sur- 
face with the annular space between the two surfaces forming a passageway 
which leads the suction flow to the compressor located at the aft end of the 
terminal fuselage taper. The turbine is located on the outside of the suc- 
tion surface where it handles the free stream boundary air produced by the 
suction process. The turbine and compressor comprise a single wheel with 
turbine blades on the outside and compressor blades on the inside. The sys- 
tem thus consists of an extremely simple and convenient arrangement of sur- 
face, ducting and turbo-machinery. 

To determine the net drag reduction of this scheme, the drag due to 
the turbo-machinery is accounted for by assuming intake and discharge pres- 
sures equal to free stream pressure. Then the drag of the compressor (suc- 
tion flow) is given by 

compressor drag = (V^ ” '^3^ 

and for the turbine 

turbine drag = (V^ - V^) 

The value of is the average velocity (based on momentum) of the boundary 
layer flow ingested by the turbine. To evaluate this parameter the asymptotic 
velocity prof i le for suction i s assumed to prevail at this station. Accord- 
ing to Reference (6) (page 231) this is given by 

V/V = 1 - exp ("C) 

CO 

where 

f = C p V y/y 

Here y denotes distance measured normal to the wall. The suction flow coef- 
ficient is taken at the optimum value = .00012 for this calculation. 

From the assumed velocity profile the average velocity is given by 
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FIGURE 30. EFFECT OF COMPRESSOR DISCHARGE MACH NUMBER ON PRESSURE 
RATIOS AND NET DRAG WITH FUSELAGE SUCTION 




0/., - X - 2 (1 - fh * 1/2 (1 - e'"’-) 

V/V “ TF 

5 ~ 1 + e ^ 


■ 25 , 


The flow in the boundary layer in terms of the suction flow is given by 


By means of these equations, the average velocity V may be found as a func- 
tion of W/W^. The results are shown In Figure (29). 

In applying these data to the turbo-machine, it is considered that the 
suction flow W is the compressor flow and the boundary layer flow W is 
the turbine flow W^. While the compressor flow is fixed by the suction re- 
quirements, the turbine flow may be chosen at will, with the choice depend- 
ing on the effect on overall performance. 


The power generated by the turbine in reducing the velocity from to 
V ; is absorbed by the compressor and utilized to reduce its drag by increas- 
ing the discharge velocity In calculating the exchange and utilization 

of power, turbine and compressor efficiencies are taken as .9 (justified by 
the low pressure ratio of the machines) and line loss recovery factors as 
•S5 (justified by the small flows involved so that ducting losses can be 
mdde very small without compromising weight or space). It is also assumed 
that all the free stream velocity head is lost in the suction process so 
that the total pressure at the compressor face Is given by x where 

is the line loss recovery factor. 

The compressor flow is fixed by the suction flow requirements, but the 
turbine flow is an independent design parameter which may be chosen for best 
overall performance. The compressor discharge velocity (or the correspond 
ing Mach number M^) is also an independent design parameter. 

The effect of is illustrated in Figure (30) which shows a typical 

variation (calculated with = W^) of the net drag (reduced skin friction 
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drag plus the momentum drags of the turbo-machine flows). It appears that 
the optimum is around .4, However, the effect on net drag is quite small 
and a value of = .2 is preferable because it entails smaller pressure 
ratios of the compressor and turbine. 

The effect of turbine flow is calculated with = .2 and with de- 
termined as a function of from the boi'r.dary layer data in Figure (29) 
with W/Wg = W^/W^. The results (Figure 31) indicate that for practical 
purposes the net drag is independent of turbine flow. Therefore, the flow 
rate may be chosen on the basis of convenience of design. The lower the 
flow rate, the higher is the required pressure ratio, two effects which are 
generally to be considered as opposing in achieving an optimum design. How- 
ever, in the present application, because flow rates and pressure ratios are 
very mild, neither quantity is of overriding importance. Therefore, it is 
considered that using a turbine flow equal to the compressor flow, leading 
to a turbine pressure ratio which is slightly higher than the compressor 
pressure ratio, is a good solution. 

Included in Figure (3l) is a line showing the drag reduction due to the 
suction only, indicating that the drag due to the turbo-machine is relatively 
small. Also shown on the figure is the result for the previous scheme of 
suction on only the lower half of the fuselage. The comparison indicates the 
greatly superior drag reduction with full fuselage suction and turbo-machine. 
Reduced drag as a function of turbo-machine efficiencies and line loss re- 
covery factor is illustrated in Figure (32). These results indicate that the 
effect of low efficiency is not significant and the effect of line loss is 
moderate. 

Of course, a fair comparison must take into account the difference in 
weight of the two schemes, the difference being favorable to the half-suction 
scheme without turbo-machine. Such a comparison involves the mission per- 
formance treated in a later section. 
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0. Wing Suction - To estimate potential drag reduction due to wing 
suction the two-dimensional results developed in Section IV-B were applied, 
on the average, to the total wing surface of the baseline configuration. 

The essential features of this scheme are indicated in Figure (33)- 

Since the suction air processed by the compressor was collected at 
variable presure, a certain ambiguity arises in defining the total pres- 
sure P|j. Accordingly, an appropriate average value is needed to treat 
compressor performance. One approach is to weign the pressure with re- 
spect to mass flow. If x represents the fraction of the flow at total 
pressure P, the average is given by 

P = / Pdx 

or if the flow is considered to consist of discrete fractions or individ- 
ual flow f i laments 

P = JiXjP. 

• I 

In order to achieve this average in practice, it is necessary to preserve 
the individual flow filaments up to the compressor blading, a procedure 
which involves transcendent ducting problems. A much more feasible ap- 
proach Is to consider the flows mixed before proceeding to the compressor, 
the ducting problem being reduced to providing a chamber where the in- 
dividual flow filaments may be mixed efficiently. In the case where this 
mixing is allowed to occur at constant area, mass and momentum conserva- 
tion imply that 

W N = ZW.N. 

where the W. represent the flow rates of the individual filaments and 

N. = N.(M.) 5 (fp) (1 + yM?) (1+^M^)^ 

fi = N(M) S (1) (1 + yM^) (1+^B^)^ 

. . -dm ■ 







In the latter expressions the M. represent the Mach numbers of the Individ- 
ual filaments at the start of the mixing process and M is the final Mach 
number of the total mixed flow. W, of course, is the total mass flow and 
is related to the W. by W = EW. . In deriving these relations it has been 
assumed that all the filaments are parallel at the start of the mixing 
process and that each have the same total temperature. 

If it is now further assumed that all of the individual filaments are 
expanded to sonic velocity at the entrance to the mixing chamber then 

N. = N = N(1) for all i 
Accordingly, the relation for N yields 

N = (N*/W) Z W. = N* 

That is, the Mach number of the mixed flow Is also sonic. It fol lows, there- 
fore, that 

W. P. A. 

^ _ _i_^ _ _j \_ 

' W PA 

where Xj 's the fraction of flow at stagnation pressure Pj > P 'S the stagna- 
tion pressure of the mixed flow, and A. and A the respective flow areas. 
Since we have assumed constant area mixing A = ZA. . Thus 



Applying this concept to the pressure distribution on the wing (Figure 2) the 
effective total pressure is determined as approximately 

p/p = .75 

OO ’ 

It is worth noting that this average pressure Is around 30 % of the basic aver 
age determined by weighing with respect to flow. 
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The momentum drag of the turbo-machine flows is determined in a manner 
similar to that described previously. The compressor provides a pressure^ 
ratio sufficient to increase the pressure of the suction flow from its ini- 
tial value P to its discharge value P„, with due allowance for the discharge 
Mach number and the line loss recovery factor. The compressor flow is fixed 
by the suction flow requirements but the turbine flow is an independent de- 
sign parameter which may be chosen at will. In making calculations, it is 
best to treat the turbine exit Mach number M 2 and the compressor exit Mach 
number M- as independent design parameters which may be chosen for best per- 
formance. The turbine entrance Mach number M^ is also a design parameter 
to some extent since the inlet of the turbine flow can be arranged to ingest 
the boundary layer flow along the fuselage or wing. 


The effect of the compressor exit Mach number M^ is illustrated in 
Figure (34) which shows curves of net drag wi th M^ = M^ = . 82 and with 
various M 2 . The value of = .2 is chosen for design since it gives very 
nearly the lowest net drag and the smal lest turbi ne flow. 

The effect of the turbine exit Mach number M 2 is illustrated in Figure 
( 35 ) which shows that a value of .3 or .4 gives the best results. The figure 
also illustrates the effect of the turbine inlet Mach number M^ , indicating 
that a serious deterioration of performance would occur with the use of 
boundary layer flow. For best overall performance (low net drag and low 
turbine flow) the turbine inlet should be free of the boundary layer and 
should handle only free stream flow. 

It should be noted that the results in Figures (34) and (35) were de- 
termined with the average total pressure at the compressor defined by 
P/P = .75 which was chosen as a realistic average value. It is interesting 
th consider the effect of this pressure on performance, illustrated by the 
results in Figure ( 36 ) calculated with M, = M^ = .82 and M 2 = .4. It is 
apparent that the average pressure has a large effect on both the net drag 

and the turbine flow. 

The effect of machine efficiencies and line loss factors are presented 
in Figure (37). The design values, shown by the circled points, result in a 
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W^, kg/s 

FIGURE 35. EFFECT OF TURBINE ENTRANCE AND EXIT MACH NUMBERS ON NET DRAG 
WITH WING SUCTION. 




FIGURE 36. EFFECT OF AVERAGE SUCTION PRESSURE ON NET DRAG WITH 
WING SUCTION 
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FIGURE 37. EFFECT OF TURBO-MACHINE EFFICIENCIES AND LINE LOSS RECOVERY 
FACTORS ON NET DRAG V/ITH WING SUCTION '' 




reduced net drag D/D^ = .928. This is less than obtained with suction on the 
fuselage for which D/D^ = .84 m. Nevertheless, the reduction is significant 
and especially when used together with fuselage suction can produce much im- 
proved mission performance as will be demonstrated in the next section. 

The physical layout of the wing-suction scheme is conceived as two sepa- 
rate symmetrical arrangements of ducting and turbo-machine, one for each of 
the port and starboard wings. The ducting of the suction flow proceeds from 
the wing tip to the wing root where the turbo-machine is located. For the 
baseline configuration the wing thickness at the root is approximately 2.1 m. 
This is large enough to accommodate the turbo-machine, the diameter of which 
is estimated to be around 1.2 m. The inlet for the turbine flow may be lo- 
cated on the bottom of the wing or on the side of the fuselage. The size of 
the inlet is very small, the area being around .2 m for each side. The 
flow exits may be located on the fuselage near the trailing edge of the wing 
root or at the trailing edge of the wing itself. This latter arrangement 
could produce some improvement in airplane performance by delaying separa- 
tion on the wing. 
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VI. MISSION PERFORMANCE 

In order to evaluate the drag reduction schemes properly it is neces- 
sary to consider the effect on the mission performance of the airplane. For 
a realistic evaluation the added weight introduced by the scheme must be in- 
cluded in the analysis. 

To keep the analysis simple, the mission is considered to consist of 
three parts. 

(1) Initial operation prior to cruise 

(2) Cruise' 

(3) Final operation after cruise 

It is assumed that range is accomplished only in the cruise part of the mis- 
sion and that all the fuel has been consumed at theend of the final opera- 
tion. 

Let us consider first the mission for the baseline configuration i.e., 
in the absence of any drag reduction. Each part of the mission is character- 
ized by a certain fuel consumption defined by the three parameters. 


Y, = F./W 
1 1 o 

where F^ , F^ and F^ are the fuel weights consumed during operations 1, 2 and 3 

and W is the gross weight of the airplane (weight at the start of operation l). 
o 

The total fuel consumed is given by 


where 


F = Y W 
o o 


= Y, . + Vj 


and the airplane weight after the mission is over is 
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W = W - F 
e o 

The range accomplished during the mission is given by 

<“'o ■ ''l> 

R = constant (L/D) log p r 

^ o 1 2^ 

where (L/D) is the maximum lift to drag ratio of the airplane. 

The constant is a function of the cruise speed and the engine perfor- 
mance, but its value is not pertinent for the present purpose since we will 
be dealing with relative mission performance. 

Consider now the airplane with a drag reduction scheme which decreases 
drag but increases weight. The decrease in drag increases the maximum L/D. 
With the assumption that the drag polar of the original airplane is parabolic 
and does not change in shape, the new maximum L/D, denoted with a prime, may 
be expressed in terms of the reduced drag D/D^ determined in the previous 
sections of this report. 

(L/D)* ^ 1 

TiT^ 2(D/DJ - 1 

The variation of this ratio with D/D^ is shown in Figure (38) If the drag 
reduction scheme did not introduce an increase in weight, this ratio would 
represent the relative increased range which would be accomplished by the 
airplane with the same fuel load. 


*tt is important to note that the effect illustrated in Figure (38) assumes 
that changes in the drag polar are due solely to changes in zero-lift drag. 
This is essentially true for the fuselage drag reduction schemes. However, 
for schemes involving wing-suction an additional influence must be accounted 
for, viz., the effect of suction on lift characteristics. For the present 
example it is estimated that an additional 2 percent increase in maximum L/D 
Is provided by the wing-suction scheme. The basis for this estimate is out- 
lined in Appendix B. 
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The weight increase of the airplane is responsible for increased fuel 
consumption required to accomplish the various parts of the mission. The fuel 
consumed in parts 1 and 3 is assumed to increase directly with the weight of 
the airplane at the start or end of these parts, denoting the case of the air- 
plane with the drag reduction by prime quantities. 


o; 


I I 

W + F 
e 

W + F 
e 


3 _ 


W 

_€ 

w 


3 e 

The drag reduction scheme introduces an increase in weight AW so that 


W = W + Alt} 
e e 


Now for completeness of the comparison two relative missions are considered, 


For the first mission the fuel load is considered to be the same or F 
In this case, there is an increase in range given by 


= F 


where 


R 


(L/D)' 


"2 ' 


(L/D) log (1 + X^) 


1 - V. 


h ' V ■ ® '''i ' V '’3' 

”2 ■ (I - V, - y'j) (I - + B) 

I 

In the latter expression for the quantity B represents the increased weight 
introduced by the drag reduction scheme, defined by 

B = AW/W 


In the second mission the range is considered to be the same or R 
In this Case, there is a reduction in the fuel load given by 


= R 
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where 


• ’ - ’'o ® (Z; (1 - Y, - Y^) 0 - Y,) - Y^) 

F/F = -p-r ''o (' - ''t' 

z‘ = (1 = - 1 

a = (L/D)/(L/D)' 

In applying these eguationSy ths weight increase is expressed by the parameter 
w which defines the weight increase per square foot of surface, 

AW = w X surface area 

In the case of complete suction on the fuselage, the surface area is taken as 

2 

the complete suction surface ir x 6.1 x 61 = 11?0 m . Iri this case the value 
of w represents the weight of the suction surface which is conceived as an 
additional surface applied over the original fuselage, the gap between the 
two surfaces forming the duct which leads the suction flow to the turbo- 
machine. The value of w does not have to be increased for extraneous ducting, 
but it must include some fraction for the structure which supports the suc- 
tion surface. The value of w may also be adjusted to account for the weight 
of the turbo-machines. However, the turbo-machine is so small that its weight 
is unimportant compared to the weight of the suction surface. 

In the case of combined fuselage suction and injection, the surface is 
taken as the suction surface on the bottom part of the fuselage or 1/2 x 1170 
= 585 m^. In this case the value of w i s larger than the weight of suction sur- 
face, since it must Include the weight required to duct and inject the Ruction 
flow on the top half of the fuselage. 

Calculations were made using an original gross weight = 3^0,000 kg 
and representative values for the fuel fractions of the original mission as 
follavs: 


= .05 
Y2 = .30 
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The values of D/D^ were taken as .8^^ for the case of complete suction with 
turbo-machine and .916 for the case of half-suction and half-injection with- 
out turbo-machine. 

The results are presented as curves of increased range (Figure 39) and 
decreased fuel load (Figure ^lO) as a function of the weight parameter. It is 
considered that the suction surface weight is around 2.5 kg/m (aluminum skin 
.076 cm thick with some allowance for structural $upport). Therefore, for the 
case of complete suction, the range increase is around l8^ or the fu'il reduc- 
tion is around lU. For the case of combined suction and injection the value 
of w is increased 20^ to account for the weight of ducting required for injec- 
tion. In this case the range increase is approximately 8% or the fuel reduc- 
tion is around 6%. 

It is apparent that complete suction offers superior mission performant?;. 
This is true even when the added weight of the turbo-machine is taken into ac- 
count. It has been estimated that the size of the wheel is less than 1.2 m in 

diameter. Since it is a single stage machine, the weight is probably around 

2 

150 kg. The added surface weight based on w = 2.5 kg/m is 2830 kg. There- 
fore, the weight of the turbo-machines increases w from 2.5 to 2.65. This in- 

crease does not reduce mission performance significantly. 

The effect of turbo-machine efficiencies and line loss recovery factor 
on mission performance (determined with w = 2.5) are shown in Figures (i+l) and 
(^2) for the case of full fuselage suction. It is evident that the effects are 
extremely mild, a direct result of the small flow involved in the suction pro- 
cess. 

Mission performance for an aircraft utilizing the wing drag reduction 
scheme whether alone or in combination with the two fuselage drag reduction 
schemes is shown In Figures {k3) and (4^). 

For the case of wing suction, mission performance is determined for the 
surface area of 900 m and a value of 0/0^ = ,928.j. (Note that the L/D factor 
in Figure 38 is increased by 2 percent to account for the improvement of lift 
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FIGURE k\. EFFECT OF TURBO-MACHINE EFFICIENCIES AND LINE 
LOSS RECOVERY FACTOR ON INCREASED RANGE WITH 
FUSELAGE SUCTION 





And 


FIGURE ^»2. EFFECT OF TURBO-MACHINE EFFIC I ENC I E:S AND LINE LOSS 
RECOVERY FACTOR ON REDUCED FUEL LOAD WITH FUSELAGE 
SUCTION 








Wing Suction 

Fuselage Half Suction and Injection 
Fuselage Suction 
Combined 1 6 2 
Combined 1 & 3 




characteristics with suction.) The appropriate value of the surface weight 
parameter is taken to be k0% higher then the basic weight 2.5 kg/m^ to ac- 
count for the ducting and turbo-machinery or w = 3.5 kg/m . As can be seen, 
the wing drag reduction scheme produces only moderate improvement, 8% in- 
crease in range or decrease in fuel load. However, when combined with 
either of the fuselage drag reduction schemes the effect is magnified. 

For the combined systems the value of D/0^ for determining the I/O 
factor is given by 

D/D^ = (0/D^) + (D/D^) - 1 

1 i 

where i = 2 or 3. The significance of the subscripts Is defined at the top 
of Figure (143). The 2 percent factor due to lift augmentation is applied 
to the combined L/D which is an approximation considered sufficiently accu- 
rate for the present purpose. 

The approximate surface areas for the combined systems to which the re- 
spective weight parameters are applied are given by 

S = S, + S. 

1 I 

The corresponding weight parameters are given by 

w S + W| S. 

„ = ^ 

There is obtained 

- 3.3; w^ = 3.0 

Using the results shown in Figures (^3) and (i»A) together with these values of 
w the performance of the various drag reduction schemes can be established. 
These are summarized in Table II below. 
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TABLE I I 

SUMMARY OF MISSION PERFORMANCE FOR 
VARIOUS DRAG REDUCTION SCHEMES 


Scheme 

w, kg/m^ 

Z Range 
1 ncrease 

% Fuel Load 
Decrease 

1 

3.5 

8 

5 

2 

3.0 

8 

6 

3 

2.5 

18 

11 

4 

3.3 

18 

11 

5 

3.0 

32 

17 


It is interesting to note that the combined performance is more than the sum 
of the individual performances. 

The effect of turbo-machine efficiencies and line loss recovery factor 
on mission performance with the combined systems is shown In Figures (i»5) and 
(46). The efficiencies and recovery factor refer to the wing suction system 
since this system is much more sensitive to these quanti ties. Furthermore, 
the recovery factor refers to the suction or compressor flow circuit since 
this is the circuit where ducting can become a problem. The results indicate 
that even with poor turbo-machine efficiencies and low recovery factor or high 
duct loss, the combined fuselage and wing drag reduction schemes produce large 
improvements in mission performance. 













VII. CONCLUDING REMARKS 


On the basis of the results presented in the previous section it is 
concluded that significant drag reduction and corresponding improvement in 
the performance of subsonic transport aircraft can be achieved by Judicious 
application of certain boundary layer control schemes. These include com- 
binations of suction and slot Injection applied to various surfaces of the 

aircraft. 

It is emphasized that the results obtained in this study are subject 
to a number of important restrictions and/or assumptions, which render the 
above cited conclusion provisional both qualitatively and quantitatively. 

These restrictions will be delineated here. 

The most fundamental assumption implicit in these studies relates to 
the possibility of preserving the laminar boundary layer state with distrib- 
uted surface suction. This assumption is particularly crucial since, if 
transition did occur under these conditions, the viscous drag would actually 
be increased above the undisturbed turbulent values. Although the feasibility 
of the technique appears to have been established in both wind tunnel and 
flight tests no applications have been made to date on either commercial or 
military aircraft (Reference l). Accordingly, our performance results have 
to be considered overly optimistic at the present time. 

On the other hand, if the aforementioned assumptions do prove to be 
valid, the performance estimates which have been, obtai ned are probably con- 
servative for a variety of reasons. As an example the wing drag reduction 
results outlined in Section V-D indicated that net reduction was quite sen- 
sitive to the average pressure at the compressor (cf: Figure 36). The per- 

formance results presented in Section VI were made wi th a single value of 
this parameter taken as P/P^ = 0.75- Higher values of this average would 
lead to substantial increases in performance. These could be achieved by 
installation of more compelx ducting possibly at the cost of increased ve- 
hicle weight. Furthermore, this parameter is a function of wing geometry. 
Thus, optimization of this parameter would be useful and could lead to im- 



proved performance. 


Additional examples of possible improvement in performance would in- 
clude application of suction on empennage surfaces. The results given in 
Appendix A are indicative of the large reduction in drag that can be achieved 
by application of this technique on airfoil sections at zero angle of attack. 

We also note that the contribution of those schemes to lift-augmentation 

has hardly been examined. Specifically, the angle of attack for which our 

wing results were made was substantially below representative values for flight 

at (L/D) . It is expected that much greater improvement in (L/D) would 

max 

result at higher angles of attack. Also, the use of slot injection on the 
wings to energize the boundary layer and prevent separation has also not been 
examined. 

In summary, the results of this study indicate a real potential in terms 
of the development of low-drag subsonic transport. However, the extent of this 
potential has not been completely established and further in-depth studies 
directed toward optimizing this potential are recommended. 
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APPENDIX A 

LAMINAR BOUNDARY LAYER WITH SIMULTANEOUS 
MASS TRANSFER AND PRESSURE GRADIENT 


The starting point for this analysis is the Von Karman momentum inte- 
gral equation for two-dimensional incompressible laminar boundary layer flow 
which takes the form (cf: Reference 6, p. 236) 


2 d0 
U T~ 
e dx 


+ (26 + 



du 
e 

^e dx 


w 


+ V 


w 


u 

e 


(A-1) 


where represents the normal velocity at the wall (negative for suction) 
and the remaining variables are defined in the conventional manner. 


Equation (A-1) can be written in the alternate form 


t' + T(H + 2) UU' = C^/2 + VU 


(A-2) 


where the transformed variables are defined by 


U 

V 

T 

H 

( ) 

X 


u /u 

e eoo 

V /u 

w e°° 

u e/v 
e“ 

6*/e 

d/dX 


goo 


. (A-3) 


Following Torda (Reference 11) we assume a velocity profile of the form 


'V 2 3 4 

u = An + Bn + Cn"^ + Dn 


(A-i») 


where 


u = u/u 


n = y/5 
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Equation (A-A) satisfies the no-slip condition u(x,0) - 0. The remaining 
coefficients are evaluated by Imposing the following additional boundary 
condi tions 


u = u^: 1^=0 ey = 6 


■■ O. ■ “ 


du .2 

e / 3 U \ 


w 


2 T 

/3iJ\ _ / 3 U \ 

( — o) = V — o) 

^ 3y^ 3y^ 

^ VI V w 


Application of (A-5) to the profile (A-4) yields 


(A- 


where 


A = (2^ + 6N + MN)/K 

B = 3(^H - 3N)/K 

C (1|M^ - 3MN)/K 
D = (3N - 6 + 2MN - 6H 

K = 18 + 6M + 

M = MR 

2 ' 

N = R U 


R = u 6/v 

goo 


3M^)/K 


also 

C^/2 = AU/R 

Substitution of these results in Equation (A-2) yields the following equation 
relating R, the Reynolds number based on boundary layer thickness, to the 
variation of mass transfer V and velocity U, 

v' = G,/UR^ + G*u" R + G,r'/R^ 



and 


= (1260K^/B^) 



= ZB^/B 
= 6B^/B 


3 

3 


{N (B^/1260K^-1) + M + B^/K} 


- (9180N^ + 7128N - 175392) + M(7290N^ - 10584N - 181728) 
+ M^(2112N^ - II388N - 86112) 

+ M^(291N^ - 3276N - 21936) + M^(19N^ - A16N - 3768) 
- M^(26N + 432) - 24M^ 

B2 = (11016N + 7128 ) + M(9126n - 6156) + M^(3114N - 8064) 

+ M^(513N - 3048) + m^(38n - 546) - 39M^ 


= (1890N^ + 26568N + 19008) + M(2898N^ + 19944N + 12672) 

+ M^(909N^ + 1368N - 19008) + M^(76N^ - 78ON - 4032) 
- M^(78N + 144) 

B^^ = 72(17N^ + 211N + 1778) + 6m(101N^ + 1007N + IO838) 

+ 2M^(38N^ + 32IN + 9108) + 3M^(11N + 936) + 216m^ 


B^ = 6N + 24 + MN 
5 

For given distributions of V and U the variation of boundary layer thickness 
can be determined from Equation (A-6) • Subsequently, all other boundary layer 
parameters (momentum thickness, skin friction, etc.) can be computed from ap- 
propriate algebraic auxiliary relations. Alternately, the condition R=0 can 
be imposed leading to 

v‘ = G^/UR^ + G^u" R (A-7) 
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Then for U specified the requisite distribution of suction needed to main- 
tain a constant boundary layer thickness can be established. 

Two features of Equations (A-6) and (A-7) are important to note here. 
First, the solution depends on the second derivative of the axial velocity 
distribution. Accordingly, in the numerical integration procedure which was 
utilized to obtain solutions a cubic-spline fit of the input velocity dis- 
tribution was employed to assure smooth variation of this parameter. 

The second feature to be noted is that the leading term on the right 
hand side of either (A-6) or (A-7) is singular at a stagnation point; i.e.: 
has no real roots. Accordingly, this method cannot be utilized to initiate 
a calculation at a stagnation point. Thus, the approximate scheme due to 
Thwaites (Reference 11) was employed for this purpose and the two methods 
matched at a small distance away from the singularity. The method of Re- 
ference (11) takes the following form. Equation (A-2) can be written 


where 


Z = P/U 


z E r 


P = 2 {c^/2U - T U (H + 2) + VT} 

Thwaites approximates the function P by 

P = 0A5 - 6T^ u' - 1,28 VT + 0.76 V^ T^ 

Real roots of this function can be found. Accordingly, the indeterminate 
form P/U at a stagnation point can be evaluated permitting integration to 

J- ■ 

be initiated there. Specifically, it can be shown that 

out 

, 2 T U P - V T P. 

2 - o 0.1 o o 2 

° u' (1 - P,) - P,/2 T„ 

o 1 o 2 o 


*This procedure represents a generalization of the method employed in the 
Karman-Pohlhausen technique to evaluate stagnation conditions in the ab- 
sence of mass transfer (see p. 210 of Reference 6). 
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FIGURE A-2. NET REDUCTION IN AVERAGE SKIN FRICTION DRAG 
AS A FUNCTION OF SUCTION PARAMETER 





where 
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1.28 + 1.52 

and the subscript o denotes values at the stagnation point. 

A computer code was developed for numerical solution of Eguations (A 6) 
and (A-7) by standard integration techniques. Some representative results 
obtained with this scheme are presented in Figure (A-l). Here the variation 
of skin friction coefficient over an airfoil surface with the indicated pres- 
sure distribution for various uniform rates of suction are presented. The 
selected pressure distribution is nominal but corresponds roughly to that 
encountered on a symmetric airfoil at a free stream Mach number of = .75 
(Reference 19). As can be seen, for suction rates moderately greater than 
the optimum value separation is suppressed with substantial reduction in skin 
friction relative to the turbulent estimates which are also shown. These re- 
sults are summarized in Figure (A-2) in terms of net reduction in average 
skin friction as a function of the suction flow coefficient. 




APPENDIX B 


LIFT AUGMENTATION DUE TO WING SUCTION 




According to Reference (9), the section lift coefficient corresponding 
to the inviscid pressure distribution shown in Figure (2) is 1.08^. Further 
it Is indicated that if viscous effects are taken into account the lift co- 
efficient is reduced to a value .779. The pressure drag coefficient is 
.00266. 




i 
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For the purpose of the present estimate it is assumed that the higher 
value of can be attained when suction is applied in accordance with the 
distribution shown in Figure (6). In the absence of suction the lower value 
prevai Is. 




It is now assumed that the two-dimensional pressure drag varies as the 

square of the lift, similar to the variation of the nominal drag polar which 

2 

L* 


has the form Cp = 


+ k Cf. In this form, the factor k is considered to be 


the sum of the three-dimensional induced drag factor and the two-dimensional 
pressure drag factor. Without improved lift, the k-factor due to pressure drag 
is given by 


Kfe 


'^1 “ =0 

P 1 


whe 


re C-. is the pressure drag coefficient and c. is the smaller of the 
Dp '-I 


stated lift coefficients. 

Assuming that the change in pressure drag is small, the. k-factor with 
improved lift is 


•^2 " 

P 2 


The improved maximum L/D is then given by 


>• 


(L/D)2 k. + 


(L/D)" *^i 4 


. If 
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where k. is the three-dimensional induced drag factor. 

I 

To evaluate the improvement, the value of k. is taken as 


k. = 


1 


1 


i ” u X Aspect Ratio ir x 7 


= .0i»5 


Substitution of the numerical data yields 


(L/D)2 ^ ^ .004i» 

(L/D)” " .0^45 + .0023 


1 .02 



